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Ph. 

fH , Abstract 

j^ , We study Schrodinger operators on an infinite quantum graph of a 

chain form which consists of identical rings connected at the touching 
points by 5-coupfings with a parameter a £ R. If the graph is "straight" , 
i.e. periodic with respect to ring shifts, its Hamiltonian has a band spec- 

CNj ■ trum with all the gaps open whenever a ^ 0. We consider a "bending" 

^ ' deformation of the chain consisting of changing one position at a single 

CTN , ring and show that it gives rise to eigenvalues in the open spectral gaps. 

We analyze dependence of these eigenvalues on the coupling a and the 
"bending angle" as well as resonances of the system coming from the 
bending. We also discuss the behaviour of the eigenvalues and resonances 

C*"- ' at the edges of the spectral bands. 

O " 
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•^ ■ 1 Introduction 
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Quantum graphs, i.e. Schrodinger operators with graph configuration spaces, 
were introduced in th e midd le of the last century [RuS53| and rediscovered three 
decades later |,GP^8l IES89 . Since then they attracted a lot of attention; they 



became both a useful tool in numerous applications and a mean which makes 
easy to study fundamental properties such as quantum chaos. We refrain from 
giving an extensive bibliography and refer to the recent proceedings volume 
|AGA| which the reader can use to check the state of art in this area. 

One of the frequent questions concerns relations between the geometry of 
a graph F and spectral properties of a Schrodinger operator supported by F. 
Put like that, the question is a bit vague and allows different interpretation. 
On one hand, we can have in mind the intrinsic geometry of F which enters 
the problem through the adjacency matrix of the graph and the lengths of 
its edges. On the other hand, quite often one thinks of F as of a subset of 
K" with the geometry inherited from the ambient space. In that case geometric 



perturbations can acquire a rather illustrative meaning and one can ask in which 
way they influence spectral properties of a quantum particle "living" on F; in 
such a context one can think of graphs with various local deformations as "bent" , 
locally "protruded" or "squeezed", etc. 

This is particularly interesting if the "unperturbed" system is explicitly solv- 
able being, for instance, an infinite periodic graph. An influence of local spectral 
perturbations mentioned above is in this setting a rich subject which deserves 
to be investigated. So far it was considered only episodically but even such a 
brief look shows that it may have properties uncommon in the usual theory of 
Schrodinger operators |KV06| . With this motivation we find it useful to start 
such a programme by discussing the influence of a "bending" deformation on a 
graph which exhibits a one-dimensional periodicity. 

To make things as simple as possible at the beginning we will not strive in 
this paper for generality and we will discuss in detail a simple nontrivial example, 
allowing for a fully explicit solution, in which the unperturbed system is a "chain 
graph" consisting of an array of rings of unit radius, cf. Fig.[Tl connected through 
their touching points. We suppose that there are no external fields. Since values 
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Figure 1: The unperturbed chain graph 



of physical constants are not important in our considerations we put h = 2m = 1 
and identify the particle Hamiltonian with the (negative) Laplacian acting as 
tpj 1-^ —■0" on each edge of the graph. It is well known that in order to get a 
self-adjoint operator one has to impose appropriate boundary conditions at the 
graph vertices. In our model we employ the so-called S- coupling characterized 
by the conditions 

n 

V',(0) = ^fc(0)=:V(0), j,keh, ^V-jW^aV^CO), (1.1) 

where n = {1, 2, . . . , n} is the index set numbering the edges emanating from 
the vertex — in our case n = 4 — and a G M U {-|-oo} is the coupling constant 
supposed to be the same at every vertex of the chain. It is important that the 
"straight" graph has spectral gap^j, thus we exclude the free boundary conditions 
(sometimes called, not quite appropriately, Kirchhoff), i.e. we assume a ^ 0. 

^A nontrivial vertex coupling is also related to the problem of approximation of quantum 
graphs by "fat graphs" of which the reader can learn more, e.g., in 13E07| or |EP08| . references 
therein, and a paper in preparation by the last mentioned authors. 



The geometric perturbation to consider is the simplest possible bending of 
such a chain obtained by a shift of one of the contact points, as sketched in 
Fig. [21 which is parametrized by the bending angle 'd characterizing the ratio 
of the two edges constituting the perturbed ring. Our aim is to show that the 
bending gives rise to eigenvalues in the gaps of the unperturbed spectrum and 
to analyze how they depend on §. At the same time the bent chain will exhibit 
resonances and we will discuss behaviour of the corresponding poles. 

The contents of the paper is the following. In the next section we analyze 
the straight chain. Using Bloch-Floquet decomposition we will show that the 
spectrum consists of infinite number of absolutely continuous spectral bands 
separated by open gaps, plus a family of infinitely degenerate eigenvalues at 
band edges. In Section[3]we will analyze the discrete spectrum due to the bend- 
ing showing that in each gap it gives rise to at most two eigenvalues. Section [4] 
describes their dependence on the bending angle as well as complex solutions 
to the spectral condition corresponding to resonances in the bent chain. In 
Section [5] we discuss further the angular dependence with attention to singular 
points where the solutions coincide with the band edges. Finally, in the con- 
cluding remarks we draw a parallel of our results with properties of quantum 
waveguides. 

2 An infinite periodic chain 

First we consider a "straight" chain Fq as sketched in Fig. [1] without loss of 
generality we may suppose that the circumference of each ring is 2tt. The state 
Hilbert space of a nonrelativistic and spinless particle living on Fp is i^(ro). We 
suppose that the particle is free, not interacting with an external potentials on 
the edges, and denote by Hq its Hamiltonian, i.e. it acts as the negative Lapla- 
cian on each graph link and its domain consists of all functions from Wi^^ (Fq) 
which satisfy the 6 boundary conditions ()l.ip at the vertices of Fq; we suppose 
that the coupling constant a is the same at each verte^o. 

In view of the periodicity of Fq, the spectrum of Hq can be computed us- 
ing Bloch-Floquet decomposition. Let us consider an elementary cell with the 
wavefunction components denoted according to the Fig. [3] and ask about the 
spectrum of the Floquet components oiHo- Since the operator acts as a negative 
second derivative, each component of the eigenfunction with energy E = k^ ^ 
is a linear combination of the functions e . The momentum k is convention- 
ally chosen positive for £^ > 0, while for E negative we put fc = Ik with k > 
(the case E = Q will be mentioned separately below). For a given E ^ Q, the 



^The coupling constant a is kept fixed and for the sake of simpHcity we will not use it to 
label the Hamiltonian neither in the straight nor in the bent case. 



wavefunction components on the elementary cell are therefore given by 



^l{x) 


^C+e^'^^ + C^e-''^", 


xei-TT/2,0] 


iPr{x) 


= C+e'^" + C^e-'^■^ 


xe [0,tt/2] 


Vl{x) 


= D+c'^'' + Dle-''''\ 


xe h7r/2,0] 


fnix) 


^D+e"'^ + D],e-"'^, 


X e [0, 7r/2] 



(2.1) 



As we have said, at the contact point the (5-coupling ()1.1|) is assumed, i.e. 

^L(0) = ^ii(0) = <^L(0)-¥'fl,(0) 

-V-lW + M^) - v'l{^) + ^fl(O) = « • ^l(O) 



(2.2) 



On the other hand, at the "free" ends of the cell the Floquet conditions are 
imposed, 

^-^(^2) = e'VL(-7r/2) M^l^) - e'Vl(-V2) 

(^fl(7r/2) = e'VL(-^/2) v'r{^/2) = c'Vl(-V2) 

with running through [— tt, tt); alternatively we may say that the quasimomen- 
tum ^9 runs through [—1/2, 1/2), the Brillouin zone of the problem. 

Substituting (|2.ip into (|2.2p and (|2.3p . one obtains after simple manipula- 
tions 

C^ • sin kir = DJ^^ • sin kir , C^^ ■ sin kir = D^ ■ sin kn , (2-4) 

where X stands for L or R, hence Cj^ = C^ and I?^ = Z?3f provided fc ^ No := 
{0, 1,2,...}. We will treat the special case fc e N later, now we will suppose 
k does not belong to N, the set of natural numbers. Furthermore, from (|2.2p 
and (|2.3p we obtain an equation for the phase factor e''', 

g2ie _ gi9 ^2 (,os fcTT + -^ sin fcTr) + 1 = , (2.5) 

V 2k J 

which has real coefficients for any k G MUiM\{0} and the discriminant equal to 

/ a \ 2 

D — [2 cos /cTT H — - sin k-K 



(^ 



2k 

We have to determine values of fc^ for which there is a G [— tt, tt) such that (|2.5p 
is satisfied, in other words, for which k^ it has, as an equation in the unknown 
e'^, at least one root of modulus one. Note that a pair of solutions of ()2.5p 
always give one when multiplied, regardless the value of fc, hence either both 
roots are complex conjugated of modulus one, or one is of modulus greater than 
one and the other has modulus smaller than one. Obviously, the latter situation 
corresponds to a positive discriminant, and the former one to the discriminant 
less or equal to zero. We summarize this discussion as follows: 



Proposition 2.1. If k"^ G ]R\{0} and k ^N, then fc^ G (t{Hq) if and only if the 

condition 

a sin fcTT 
cos KTT H — • — - 
4 h 

is satisfied. 



< 1 (2.6) 



In particular, the negative spectrum is obtained by putting fc = iK for k > 
and rewriting the inequahty (|2.6p in terms of this variable. Note that since 
sinhx 7^ for all x > 0, it never occurs that sinfcTr = for k G iM+, the 
positive imaginary axis, thus there is no need to treat this case separately like 
for fc e R+, cf. (Hill) above. 



Corollary 2.2. If k > 0, then —k? e (T{Ha) if and only if 

<1. (2.7) 



a sinh ktt 

cosh KTT H • 



4 K 

Let us finally mention the case fc G N left out above. It is straightforward 
to check that k"^ is then an eigenvalue, and moreover, that it has an infinite 
multiplicity. One can construct an eigenfunction which is supported by a single 
circle, which is given by '4^{x) = svcikx with x G [0,7r] on the upper semicircle 
and ^p{x) — — sin kx with x G [0, tt] on the lower one. 

Remark 2.3. The condition (|2.6|) reminds us of the corresponding condition in 
the Kronig-Penney model with the distance between the interaction sites equal 
to TT, cf. [AGHHJ ■ the only difference being that the coupling constant is halved, 
ia instead of a. In contrast to that, the point spectrum of the KP model is 
empty. These facts are easy to understand if we realize that our model has 
the up-down mirror symmetry, and thus H^ decomposes into a symmetric and 
antisymmetric part. The former is unitarily equivalent to the KP model with 
modified coupling, the latter corresponds to functions vanishing at the vertices, 
having thus a pure point spectrum. Looking ahead, we remark that the bending 
perturbation breaks this mirror symmetry. 

Finally, in the case E — Q we get in the similar way the equation 

e2i«_ei«(2 + ^)+l = 0, (2.8) 

replacing (|2.5p . whence we infer that G (t{Hq) if and only if 

<1, (2.9) 



an 
1 + - 



hence zero can belong to the continuous part of the spectrum only and it happens 
iff a £ [— S/tt, 0]. In conclusion, we can make the following claim about a{Ho). 

Theorem 2.4. The spectrum of Hq consists of infinitely degenerate eigenval- 
ues equal to n^ with n G N, and absolutely continuous spectral hands with the 
following properties: 

If a > 0, then every spectral band is contained in an interval (n^, {n + 1)^] with 
n £ No :— N U {0}, and its upper edge coincides with the value (n + 1)^. 

If a < 0, then in each interval [n^, (n + 1)^) with n G N there is exactly one 
spectral band the lower edge of which coincides with n^ . In addition, there is a 



spectral band with the lower edge (being the overall spectral threshold) equal to 
—K^, where k is the largest solution of 



a sinh ktt 

cosh KTT H 



4 K 



1. (2.10) 



The position of the upper edge of this band depends on a. If — S/tt < a < 0, 
then it is equal to k^ where k is the solution of 

a sinfcTT 

cos k-K + — ■ ; = — 1 

4 k 

contained in (0, 1). On the other hand, for a < — S/tt the upper edge is negative, 
~k} with K being the smallest solution of (|2.10p . and for a — —S/tt it equals 
zero. 

Finally, a{Ho) = [0, +oo) holds ifa — 0. 

Proof. The degenerate bands, in other words, the eigenvalues of infinite multi- 
plicity, were found already and it is straightforward to check that no other eigen- 
values exist. The continuous spectrum can be in view of Remark 12. 31 treated as 
in AGHHJ , nevertheless, we sketch the argument not only to make the paper 
self-contained, but also in view of next sections where some ideas and formula 
of the present proof will be used again. 

Consider first the positive part of the continuous spectrum. The condi- 
tion (|2.6|) clearly determines bands with one endpoint at n^, n G N, where the 
sign of a decides whether it is the upper or lower one. If a < 0, the pres- 
ence of a band in (0,1) depends on |a|. Denoting g{x) :— cosxtt -|- j • Sii^i 
we want to show that B := {x € (0,1) : \g{x)\ < 1} is either empty or an 
interval with zero as its edge. It is obvious that g maps (0, 1) continuously into 
(—00,1); we will check that g{xo) = —1 implies g'{xo) < 0. We notice first 
that the premise implies cosa^o''' = ^1 ~ J ' ^"^^°^ ; taking the square of this 

relation we find after simple manipulations that sinxovr = —2 f -^ + ^^ j 

and cosxoTT = (j^ ^ j (^ + ^ j . Evaluating g'{xo) and substituting 

these expressions we get 

a-K ( sin7rao\ 
Axq \ 7ra;o / 

These properties together with the continuity of g imply that if B is not empty, 
then it is an interval with the left endpoint zero. It is also clear that B is non- 
empty iff g(0+) > — 1 which gives the condition a > —S/tt. On the contrary, B 
is empty if a < —S/tt and the borderline case a = —S/tt was mentioned above. 
Let us next focus on the negative part using g{x) :— coshxTr-l- j ■ '^'"^^^ and 
ask about _B := {x e (0, oo) : \g{x)\ < 1}. It is easy to check that g{x) = —1 
iff tanh ^ = 45- and g{x) — 1 iff coth ^ — ^- It implies that there is exactly 
one xi such that g{xi) — 1, and that the equation g(x) = — 1 has one solution 



X-i in the case a < — S/tt and no solution in the case a S [— S/tt, 0). Since 
obviously < x_i < xi and ^(0+) := lima;^o+ dix) — 1 + a7r/4, we infer that B 
is a bounded interval. Its closure contains zero iff a £ [~8/7r, 0) because then 
5(0+) G [—1,1). In such a case the lowest spectral band is the closure of i? U B, 
otherwise it is the closure of B only. D 

3 The perturbed system 

3.1 General considerations 

Let us suppose now that the straight chain of the previous section suffers a 
bending perturbation as shown in Figure [21 We call the perturbed graph F^; 
it differs from Fo by replacing the arc lengths tt of a fixed ring, conventionally 
numbered as zero, by tt ± i9. The bending angle § is supposed to take values 
from (0, tt), regardless of the fact that for d > 27r/3 it is not possible to consider 
F^ as embedded in the plane as sketched — one can certainly realize such a 
"bending" in an alternative way, for instance, by deforming the selected ring. 

The state Hilbert space of the perturbed system is i^(F^) and the Hamil- 
tonian is H{) obtained by a natural modification of i/p; our aim is to determine 
its spectrum. Since F^ has the mirror symmetry w.r.t. the axis of the zeroth 
ring passing through the points x = ^ (tt ± ^) , the operator H^} can be reduced 
by parity subspaces into a direct sum of an even part, H^ , and odd one, H^\ 
for the sake of simplicity we drop mostly the subscript -d in the following. 

All the components of the wavefunction at energy fc^ are linear combinations 
of e^'*"'^ . As we have said we use the ring labelling with zero corresponding to 
the perturbed one; the mirror symmetry allows us to study a half of the system 
only, say, with non-negative indices. The wavefunction on each ring will be a 
pair of functions ipj and Lpj , where j is the circle index, ij^j corresponds to the 
upper semicircle and Lpj to the lower one, 

^,{x)^C+c'^'^+C-e-''^\ xe[0,7r], 
^^{x) = i?+e'^-" + DJe-''^^ x e [0, tt] 

for j e N. The situation is different in the case j = 0, where the variables run 
over modified intervals, 

V'o(x)=CoV" + Co-e-''=^ xe 

(3.2) 
Mx)^D+e"'^ + D^e-"'^, xe 

There are (5-couplings with the parameter a in the points of contact, i.e. 

^,(0) = ^,(0) V,(7r)=¥',W (3.3) 
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and 



V',(0) = ^,_i(7r) (3.4) 

^'^{0) + ^;(0) - ^/y^_,{n) ~ ^;_i(7r) = a • V,(0) (3.5) 



Substituting (|3.ip into (|3.3p we obtain 



C^ ■ sin fcTT = D^ ■ sin fevr and C • sin /ctt = D- ■ sin /ctt , 

thus for fc ^ No we have C^ — D^ and CJ^ — DJ . The case fc G No can be 
treated in analogy analogously with the "straight" case: it is easy to see that 
squares of integers are infinitely degenerate eigenvalues and the eigenfunctions 
can be supported by any ring, now with the exception of the zeroth one. From 
now on, we suppose fc ^ No. 

Using the coupling conditions (|3.4[) and (|3.5[) . we arrive at a "transfer ma- 
trix" relation between cocfhcients of the neighbouring rings, 






(3.6) 



M 



valid for all j > 2, which yields 



^^ ^ - AP-' ■(% ) . (3.7) 



CJ / V c. 



It is clear that the asymptotical behavior of the norms of {C^,C^ ) is de- 
termined by spectral properties of the matrix M. Specifically, let {C^,C^)'^ 
be an eigenvector of M corresponding to an eigenvalue /i, then |/x| < 1 
(|/i| > 1, |/i| = 1) means that ||(Cj^, C~)'^|| decays exponentially with respect 
to j (respectively, it is exponentially growing, or independent of j). 

The wavefunction components on the j-th ring for both H (as well as on 
the (-j)-th by the mirror symmetry) are determined by Cf and C^ , and thus 
by {Ci , Cf)^ by virtue of ()3.7|) . If (Cj^, Cf )-^ has a non- vanishing component 
related to an eigenvalue of M of modulus larger than one, it determines neither 
an eigenfunction nor a generalized eigenfunction of H . On the other hand, 
if (C^,Cj")-^ is an eigenvector, or a hnear combination of eigenvectors, of the 
matrix M with modulus less than one (respectively, equal to one), then the 
coefficients C- determine an eigenfunction (respectively, a generalized eigen- 
function) and the corresponding energy E belongs to the point (respectively, 
continuous) spectrum of the operator H . To perform the spectral analysis of 
M, we employ its characteristic polynomial at energy k'^, 

X^ -2x(coskTT+ —smkn) +1, (3.8) 

V 4fc / 



which we have encountered already in the relation (|2.5[) : it shows that M has 
an eigenvalue of modulus less than one iff the discriminant of (13. 8p is positive, 



I.e. 

a 



cos fcTT H sin kiT 

Ak 



>1, 



and a pair of complex conjugated eigenvalues of modulus one iff the above 
quantity is < 1. In the former case the eigenvalues of M are given by 



— - sin fcTT ± -\ / I c 
4fc V V 



A 1 = cos k-K + -rr sin kn ± \l { cos kir + — sin fcyr J — 1 , 

satisfying A2 — Aj^^, hence A2 < 1 holds if coskn + -^ sinfcTr > 1 and Ai < 1 if 
this quantity is < —1. Moreover, the corresponding eigenvectors of M are 

/ a p— ifc-TT 

V ^1,2- (.l + ilfcje 

Remark 3.1. Comparing to (|2.6p we see that the perturbation does not affect 
the spectral bands, and also, that new eigenvalues coming from the perturbation 
can appear only in the gaps. These facts are obvious, of course, from general 
principles. Using the natural identification of L^(ro) and L^(r^) we see that Hq 
and H^ differ by a shift of the point where a boundary condition is applied, hence 
their resolvent difference has a finite rank (in fact, rank two). Consequently, 
their essential spectra coincide and each spectral gap of Hq contains at most 
two eigenvalues of H^, see jWel Sec. 8.3, Cor. 1]. 

3.2 Spectrum of H+ 

The operator H^ corresponds to the wave functions even w.r.t. the symmetry 
axis, hence we may consider a half of the graph with the Neumann conditions 
at the boundary (i.e., the points A, B in Figure [21), 

At the contact point of the zeroth and the first ring (denoted by C) there is a 
(5-coupling with the parameter a, 

^o(^) = '^o(7r)=^i(0) (3.9) 

^[{0) + ifi'M ~ V^(7r) - ^[,(7r) = a ■ M^) (3-10) 

Substituting to these conditions from (|3.ip and (|3.2p and using the equality 
(p'i{0) = V'i(O), we obtain (C^, C^)^ up to a multiplicative constant, 

_l_ V / cos fcTT+cos kit I j f T _ a(cosfc7r+cosfc)?) 

O"!^ \ / sink-K \ 2fc sin fcTT 



Ct J \ cos fcTT+cos fci? _ j ( 1 _ a(cosfc7r+cosfci?) 

sin k-K \ 2k sin hrr 

The right-hand side is well defined except for sin kir = 0, but this case has been 
already excluded from our considerations; we know that for fc G N the number 
/c^ is an eigenvalue of infinite multiplicity. 



Following the above discussion k"^ G ap{H^) requires that the vector (C^, Cf )-^ 
is an eigenvector of M corresponding to the eigenvalue A of the modulus less 
than one. Using the above explicit form of the eigenvectors and solving the 
equation 



cos fcTT+cos fci? 4- i I 1 — a(cosfc7r+cosfcT?) \ _Q_ — ifcir 

sin kn I 2k sin kTZ J 4ik 

cos fcTT+cos fci? _ j ( 1 _ a(cosfc7r+cosfci?) \ \ _ /'l j ^\ ^ikTr 

sin /cTT \ 2fcsinfc7r 



A-(l + 4)e" 



= 



we arrive at the condition 



(cos k-d + cos kir) ■ — sin kn ± \ ( cos /stt + — sin fcyr j — 1 = sin'^ /ctt , 

with the sign given by the sign of cos kir+jj: sin kir. Since sin kir ^ 0, the second 
factor at the Ihs is also nonzero and the last equation is equivalent to 

7 a 7 sin^ /cTT /„ 1 1 N 

cos fcu = — cos KTT H ; (3.11) 

^ sin fcTT ± -W (cos fcTT + ^ sin fcyr) — 1 

for the sake of brevity we denote the expression at the rhs by /(fc). 

The relation (j3.1ip is our main tool to analyze the discrete spectrum and 
we are going to discuss now its solutions. We start with an auxiliary result 
noting that, as a consequence of Theorem 12.41 the set of positive k for which 
the inequality | cosfcyr + ^ sin/cTrj > 1 is satisfied is an infinite disjoint union 
of closed intervals. We denote them /„ with n e N and recall that n G J„. If 
a > we denote by /q the interval with the edge at zero corresponding to the 
non- negative part of the lowest spectral gap of Hq . 

Proposition 3.2. The function f introduced above maps each In\{n} into the 
interval (-1,1) U {(-1)"}. Moreover, f{x) = (-1)" holds for x G /„\{n} iff 
|cosa;7r + ^sinx7r| = 1, and limj.g/„_j.^„ /(x) = (-1)"+^ 

Proof. According to (j3.1ip . the function / is continuous in each interval /„\{n}, 
thus it maps the interval 7„\{n} again to an interval. The claim then follows 
from the following easy observations. First, f{x) = (—1)" iff x is the non- 
integer boundary point of J„ (if a < and |a| is sufficiently large, the left edge 
of Ji is moved to zero and one checks that lima;_jo /(a;) = —1). Furthermore, 
for all X G In\{n} we have f{x) ^ (—1)""^, and finally, linix^n.xein /(•^) = 
(-1)"-!. ' " D 

Proposition 13.21 guarantees the existence of at least one solution of (|3.1ip in 
each interval /„\{n}, except for the case when -d satisfies cosn^ — (—1)"^^, or 
equivalently, except for the angles d = "^^~^^ 7r, t = 1, . . . , [^^^] . Later we 
will show that for these angles there is indeed no solution of the equation (|3.1ip 
in /„\{n}, while for the other angles in (0,7r) there is exactly one. 
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In a similar way one can proceed with the negative part of the spectrum. If 
k = IK where k > 0, the condition (|3.11[) acquires the form 

1 n 1 sinh KTT ,„ , „, 

cosh KW = — cosh KTT — _ (3.12) 

^ sinh KTT ± y (cosh nn + -^ sinh ktt) — 1 

where the upper sign in the denominator refers to cosh kit + ^ sinh ktt > 1, and 
the lower one to coshKTr + -^ sinhKTr < —1. Let us denote the rhs of p.l2p by 
/(k), then we have the following counterpart to Proposition 13.21 

Proposition 3.3. If a > 0, then /(k) < — coshKi? holds for all k > and 

"d G (0, n). On the other hand, for a < we have 

// limK_,o fcoshK7r + ^sinh^Tr) < —1, then there is a right neighbourhood of 
zero where f{x) = —1 — C{a)x'^ + o{x'^) with the constant explicitly given by 

C{a) := f^+{^ + ^[9fY ^s^yj J 7^2. Moreover, /(k) = -1 holds for 

K > iff cosh KTT + ^ sinh ktt = — 1 . 

The interval {k : coshKTr + ^ sinh^TT > 1 A k • tanhKTr < — a/2| is mapped by 
the function f onto [1, +oo). 

If K tanhKvr > —a/ 2, then /(k) < — cosh^i^ holds for all k > and i9 G (0, tt). 

Proof. The statement for a > is obvious, assume further that a < 0. The first 
claim follows from the Taylor expansions of the functions involved in /, the last 
uses the equality cosh k — sinh k = 1. The set determined by the conditions 
cosh KTT + -^ sinh KTT > 1 and k ■ tanhKTr < — q;/2 is obviously an interval and 
/ is continuous on it. Since coshKTr + ^sinh^Tr ~ 1 implies /(k) = 1 and 
for Kg • tanhKoTT = —a/2 it holds hm^^^- / = +oo, the second claim follows 
immediately. Finally, if k ■ tanhKTr > —a/2, then coshKTr + ^sinhKTr > 1 

and ^ sinh KTT ± y (coshK7r+ ^sinhKTr) — 1 > 0, thus /(k) < — coshKTr < 
— coshKi? holds for all k > and i? G (0, tt). D 

In particular, the first claim concerning a < together with the continuity 
of / implies that if the set { k ; cosh ktt + ^ sinh ktt > 1 } is nonempty (and 
thus an interval), the graph of / on this set lies below the value -1 touching it 
exactly at the endpoints of this interval. 

Corollary 3.4. If a > 0, then iJ+ has no negative eigenvalues. On the other 
hand, for a < the operator H^ has at least one negative eigenvalue which 
lies under the lowest spectral band and above the number —Kq, where kq is the 
(unique) solution of k ■ tanhKvr — ~a/2. 

Proof. The eigenvalues are squares of solutions to the equation coshKi? — f{K). 
The absence of negative eigenvalues for a > follows directly from the first 
claim in Proposition 13.31 The same proposition implies that there is exactly 
one interval mapped by / onto [1, +oo), hence there is at least one solution of 
coshKi? = f{i^) in this interval. D 
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3.3 Spectrum of H and a summary 

The operator H~ which corresponds to the odd part of the wavefunction can 
be treated in an analogous way. The boundary conditions on the zero circle are 
now Dirichlet ones, 



One can easily find the spectral condition, 

sin fcTT 
— cos fcty = — cos KTT H ; (3.f3) 

Jl sin fcTT ± W (cos fcTT + ^ sin fcyr) — 1 

in comparison with (|3.11[) corresponding to H^ there is a difference in the sign 
of the cosine on the left-hand side. Since we already know the behaviour of the 
right-hand side, cf. Proposition l3.21 we can infer, similarly as for 7?+, that there 
is at least one solution of (j3.13|) in each interval J„ except for the case when 
-cosn'd= (-1)"~\ i.e. when i9 = ^^tt, £= 1,..., [f]. 

Following the analogy with the symmetric case further we can employ Propo- 
sition 13.21 to conclude that in each interval /„ there is at least one solution of 

— cosk-d — /(k). The only exception is the interval /i for a < 0: for \a\ suffi- 
ciently small it holds — cosk-d < f{k) in the whole /i; we will comment on this 
situation in more detail in the next section devoted to resonances. The negative 
part of the point spectrum of H^ is determined by the condition 

1 n 1 sinh KTT /„ , -N 

— cosh K?7 = — cosh KTT — , (3.14) 

-^ sinh KTT zt \/ (cosh kit + -^ sinh ktt) — 1 

where we set k = in for k S M"*". It follows from Proposition 13.31 that ()3.14|) has 
a solution for negative a only, and it happens if (i) the positive spectral gap 
touching zero extends to negative values, and (ii) the bending angle i!) is small 
enough. In other words, if there is a number kq solving cosh ktt + ^ sinh ktt = 

— 1, the energy plot w.r.t. ^ obtained as the the implicit solutions of (|3.14p is a 
curve departing from {d, E) = (0, —n^)\ in the next section we will show that it 
is analytic and following it one arrives at the point (i?, E) = (tt, 1). 

Let us summarize the discussion of the discrete spectrum. We have demon- 
strated that for each of the operators H^ there generally arises at least one 
eigenvalue in every spectral gap closure. We have also explained that such an 
eigenvalue can lapse into a band edge equal to n^, n S N, and thus be in fact 
absent. The eigenvalues of iJ"*" and H~ may also coincide, in this case they 
become a single eigenvalue of multiplicity two. One can check directly that it 

happens only if 

, , a 

k ■ tan kn = — . 
2 

The study of the resonances, performed in the next section, will help us to 

find more precise results concerning the number of eigenvalues. We will show 
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that there are at most two of them in each spectral gap. However, to make the 
explanation clearer, we refer already at this moment to the Figs.|4H6]illustrating 
the numerical solution of the spectral condition for different signs of the coupling 
constant, as well as the resonances of the system. 

4 Resonances and analyticity 

Proceeding further with the discussion we want to learn more about the angle 
dependence of the perturbation effects. First we note, however, that the added 
eigenvalues are not the only consequence of the chain bending. One has to 
investigate all solutions of (|4.1[) . not only the real ones which correspond to 
ap{H'^), but also complex solutions describing resonanceSL] of H . 

Proposition 4.1. Given a non-integer k > 0, the conditions (I3.1ip and (|3.12p 

for H , respectively, are equivalent to 

— (l±cosfc^cosfc7r)(±cosfc^ + cosfc7r) = sin/c7r- (l±2cosfc^cosA:7r + cos^ kd) . 

(4.1) 

Proof. First we note that changing the square root sign in denominator of (j3.1ip 
does not give rise to a real solution. Indeed, if the sign of the right-hand side 
of p. lip is changed, the obtained expression is of modulus greater than one, 
hence it cannot be equal to cos k'd. This further implies that one need not 
specify the sign in the denominator of p. lip by the sign of cos kir + -^ sin /ctt, 
and therefore we can express the square root and subsequently square both sides 
of the obtained relation. After simple manipulations, we arrive at (14. ip : note 
that for all k G ]R+\N, the denominator of p. lip is nonzero. The equivalence 
of p. lip and (|4.ip for k £ C \ N is obvious for (|4.ip considered with the complex 
square root, i.e. without restrictions on the sign in the denominator. The 
argument for H^ is analogous. D 

Now we are ready to state and prove the analyticity properties. Since the 
cases of different symmetries are almost the same, apart of the position of the 
points where the analyticity fails, we will mention the operator H^ only. 

Proposition 4.2. Curves given by the implicit equation (|4.ip for H^ are 

analytic everywhere except at (f?, fc) = ( ""'""'""^^ Tr, n), where n G N, ^ G No, 
^ < [^^^Y^l . Moreover, the real solution in the n-th spectral gap is given by a 
function iS i-^ k which is analytic, except at the points "^^^ tt. 

Proof. First we will demonstrate the analyticity of the curves i!) i-^ k £ C This 
is easier done using equation p.lip : we have to prove that at each point {d, k) 
solving the equation G(i?, fc) = with 

G(i9, k) — ~ cos k'd — cos kn -\ -j^^^^^^^=^^^^^^^ 

■^ sin kiT ± y (cos fcyr + ^ sin fcTr) — 1 



•^The notion of resonance in the chain-graph system can be introduced in different, mutually 
equivalent, ways similarly as in IEL07| . 
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the derivative ^ is nonzero. We have ^ — k ■ sin(fcz?) — iff sink-d — 0, 
i.e. k-d — rmr, m G Z. This implies G{{),k) — (—1)™+^ — cosfcTr, and since 
G{'0, fc) = should be satisfied, k is an integer of the same parity as to + 1. For 
A: e N, G is not defined and we use (|4.ip : it is easy to check that any solution 
{d, k) of (|4.ip with fc G N corresponds to 

„ k + l-2(. „ ^r . \k + l 

'& = —— TT, £eN,£< 

k 

To prove that real solutions are analytic functions, it suffices to check that, 
except at the points {(p,k) — ( "+1-2^ ^^ ^^'^^ fgp each (??, fc) solving F{d,k) = 
with 

F^-djk) := a(l + cos kiD cos kir) (cos kd + cos kir) 
—2kainkTT ■ (1 + 2 cos fci9 cos /ctt + cos^ k-d) 

it holds 4jr 7^ 0. Computing the derivative one obtains an expression which can 
be cast into the form 

2 sin^ kiT ■ {1 + 2 cos kd cos kir + cos^ fci?)^+ 

+ a ■ [7r(cos k-d + cos kn)'^ + sin^ kn ■ (cos /ci? + cos kn)'^ + 
+■& sin^ kn sin^ /^(^(l + cos k{'K — 1})) + {t: — ■&) sin^ fevr sin^ /c'i9(l + cos kd cos fevr)] 

This is always non-negative, and vanishes iff 

(cosfcTT = 1 A coski} = —1) V (cos/cTT = — 1 A cosk-d = 1) . 

i.e. iffkEZ and k-n = kd+ (2£ — l)7r, i' e Z, proving this the sought claim. D 

The resonance dependence on the bending angle ■& is again visualized on 
Figs. |4H6] where the real parts are shown; the imaginary parts corresponding to 
the situation of Fig. 2] are plotted on Fig. [7] below. 

5 More on the angle dependence 

The above results raise naturally the question about the behaviour of the curves 
at the singular points [3, k] = [ "+^~^V , n] with neN,£ eN,£ < [2±i] , where 
they touch the band edges and where the eigenvalues and resonances may cross. 
Now we are going to examine the asymptotic expansion at these points and to 
look how many curves "stem" from them. 

Consider again first the part _ff+. Let fco G N and ^0 '■— "+-'-~^^ 7r for some 
£ eN, and put 

fc:=fco + e, ■d:=i}Q + S. 

After substituting into (|4.ip with the plus signs and employing Taylor expan- 
sions of the cos and sin functions we arrive at the relation 

I {k^S^ + 4fc3^o'5'£ + 6ki^lS\^) - fco^^e' = 0{Se^) + 0{e^) + 0{Sh) . 
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Using the theory of algebroidal functions and Newton polygon, we find that in 
the neighbourhood of (fJoj ^o), the asymptotical behaviour of solutions is given 
by the terms of the order S^ and e^. In other words, up to higher-order term 
we have f fcg^^ = fcoTr^e'^, and therefore 

Note that a G R, feg > 0, (5 G M, i.e. only e may be complex here, hence the last 
equation admits exactly three types of solutions: 

• e ~ \/j^S^ (a- I'eal solution corresponding to the spectrum) 

• e = e^'3^ y/jit'^^^'^ (imaginary solutions corresponding to resonances) 

Let us remark that since (j4.ip has a symmetry with respect to the complex 
conjugation of fc, the imaginary solution come in pairs. This is why we find 
pairs of curves outside from the real plane, conventionally just one of them is 
associated with a resonance. 

Returning to properties of eigenvalues in a fixed spectral gap, we have so far 
demonstrated that each real curve describing a solution of (|4.ip is a graph of 
a function analytic except at the singular points, cf. Proposition 14.21 Further- 
more, at each singular point only one pair of branches meets (with respect to 
the variable d); it follows that there is exactly one solution in each spectral gap 
closure. Assuming for definiteness a > we can say that the complete graph of 
solutions of (|4.ip has the following structure: 

• It consists of curves that are analytic and not intersecting, except at the 
points (i?,fc) = (S+Ip2^7r,n), where n e N, i £ N, i < [^]; these are 
the only ramification points. 

• The real curves branches join the points ( ""'""'""^^ Tr, n) and ( "+i-2£-2 ^^ ^-^^ 
i.e. the consecutive points on the lines k = n Q N. 

• The curves branches outside the plane $5(fc) = join the points {—zj'^, n — £) 
a^d ( ^_"^_j^ 7r, n—£~l), i.e. the consecutive points laying on the hyperbolas 
(^ + ^) ■ fc == n • TT, fc e R, n e N, n odd, cf. Fig. El 

Furthermore, we have seen that the behaviour of eigenvalues in vicinity of the 
singular points is the following, 

V 4 TT 

and this is valid for in the particular case iSq — 0, ko £ N, as well provided the 
band edge fco is odd. 

However, i?+ has an eigenvalue near i?o = also in the gaps adjacent to 
even numbers. In these cases the curve starts at the point (0, kg) for kg being 
the solution of | cosfevr -I- ^ sin/cTrj = 1 in (n,ri + 1), n even. The asymptotic 
behaviour of k for •& close to zero is then different, namely: 
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Theorem 5.1. Suppose that n E N is even and fcg is as described above, i.e. 
kg is the right endpoint of the spectral gap adjacent to n^ . Then the behaviour 
of the solution of (|4.ip in the neighbourhood of{0,ko) is given by 

fc = fco - Cfc„,„ • 7?4 + 0(^5) ^ 

where Cko,a := ^ ■ (f ) (koT^ + sinfcoTr)" • 

Proof. The argument is straightforward, it suffices to use Taylor expansions 

in ^i^. n 

The analogous asymptotic behaviour applies to fc^, the energy distance of 
the eigenvalue from the band edge is again proportional to i}* in the leading 
order. Notice that this is true in any spectral gap, but of course, the error term 
depends in general on the gap index. 

We refrain from discussing in detail the odd part H^ of the Hamiltonian. 
The corresponding results are practically the same, the only difference is that 
the roles of the even and odd gaps are interchanged. 

Most of what we have discussed above modifies easily to the case of attractive 
coupling with the obvious changes: for a < the spectral gaps lay now below 
the numbers n'^, n € N. Of particular interest is the spectral gap adjacent to 
the value one, because with the increase of |a| its lower edge moves towards 
zero and may become negative for \a\ large enough. The even part H~^ has 
similar properties as before: the eigenvalue curve goes from (0,1) to (7r,fco), 
where fco S (0, 1), and there two complex conjugated branches with ^{k) > 
one of which describes a resonance. 

However, the odd part H~ requires a more detailed examination. We know 
that there is an eigenvalue curve going to the point [n, 1]. If the entire spectral 
gap is above zero, this curve joins it with [0, fcg], where k^ is the lower edge of 
the gap. On the other hand, if \a\ is large enough the eigenvalue curve starts 
from [0, —kq], where —Kg is again the lower gap edge; to show that even in this 
case the curve joins the points [0, —kq] and [tt, 1] analytically, it suffices to prove 
that the solutions of (|4.ip with the negative sign preserves analyticity when it 
crosses the line k"^ = 0. 

The spectral condition (|3.13p for H^ is valid for k ^ 0. If we put all terms 
to the left-hand side denoting it as Q^{{},k), i.e. 

Q~ I'd, k) = — cos kd + cos kir ■ 



f^ sin kn ± y (cos kn + -^ sin fc-Tr) — 1 



with the sign in the denominator properly chosen, we have lim/j^o Q~ {-d, k) k~^ = 
for ^ = 0,1 while for I = 2 the limit is real- valued and non- vanishing. It follows 
that to find the behaviour at the crossing point one has to examine the function 
given implicitly by G{'d, k) ^ 0, where 



S~(i?,fc) 



G(^,fe) 



for fc ^ 



g-(s,k) 



limfc^o — j^ for fc = 
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This is continuous and it can be easily checked that it has continuous partial 
derivatives with respect to z? and k in the neighbourhood of any solution of 
G{'d, k) = with fc = 0. In particular, the derivative w.r.t. d equals k~^ sin fci? 
for all fc ^ 0, thus at a point [i^q, 0] solving €!{•&, fc) = we have 

dGido,0) sinfc^o 

av k^o k 

in other words, the solution of G'(^, fc) = is analytic also at the point [^Oj 0]. 
Needless to say, this claim which he have checked directly here can be obtained 
also by means of the analytic perturbation theory Ka66] . 

Finally, note that by Proposition 14.21 the solutions of (|4.1D with both the 
positive and negative signs are analytic in the whole open halfplane 5R(fc) < 0, 
and consequently, no resonances curves can be found there. 

6 Concluding remarks 

We have reasons to believe that the spectral and resonance properties due geo- 
metric perturbations of the considered type hold much more generally. In this 
paper we have decided, however, to treat the present simple example because it 
allowed us to find a rather explicit solution of the problem. 

The problem can be viewed from different perspectives. As an alternative 
one may interpret the chain graph as a decoration of a simple array-type graph, 
or if you wish, the Kronig-Penney model, in the sense of [AIOO] and [Ku05| . 
The results of the paper then say that a local modification of the decoration can 
produce a discrete spectrum in the gaps and the other effects discussed here. 

It is also interesting to draw a parallel between the quantum graphs discussed 
here and quantum waveguides, i.e. Laplacians in tubular domains. Although 
the nature of the the two system is very different, they nevertheless share some 
properties, in particular, the existence of bound states below the essential spec- 
trum threshold due to a local bend. This effect is well studied for Dirichlet 
quantum waveguides where it is known for a gentle bend the binding energy is 
proportional to the fourth power of the bending angle |DE95) . i.e. it has exactly 
the same behaviour as described by Theorem 15. II 

Bent quantum waveguides with mixed (or Robin) boundary conditions were 
also studied |Ji06| and it was shown that the effect of binding through bending is 
present for any repulsive boundary. In our case an eigenvalue below the lowest 
band exists whenever a 7^ which inspires another look at the waveguide case. 
It appears that the argument of |Ji06| works again and proves the existence 
of curvature- induced bound states in all cases except the Neumann boundary 
which is an analogue of the case a = here. 
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Figure 2: A bent graph 



20 




Figure 3: Elementary cell of the periodic system 
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Figure 4: The spectrum of H~^ as a function of -d for repulsive coupling, a = 3. 
The shaded regions are spectral bands, the dashed lines show real parts of the 
resonance pole positions discussed in Sec. [H 
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Figure 5: The spectrum of H in the same setting as in Fig. |4] 
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Figure 6: The spectrum of iJ as a function of i? for attractive coupling, a — —3. 
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Figure 7: The imaginary parts of resonance pole positions in the same setting 
as in the previous picture; for the sake of lucidity only the curves corresponding 
to _ff+ are plotted. 
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